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We calculate the charge and spin Drude weight of the one-dimensional extended Hubbard model 
with on-site repulsion U and nearest-neighbor repulsion V at quarter filling using the density-matrix 
renormalization group method combined with a variational principle. Our numerical results for the 
Hubbard model {V = 0) agree with exact results obtained from the Bethe ansatz solution. We obtain 
the contour map for both Drude weights in the f/y-parameter space for repulsive interactions. 
We find that the charge Drude weight is discontinuous across the Kosterlitz-Thouless transition 
between the Luttinger liquid and the charge-density-wave insulator, while the spin Drude weight 
varies smoothly and remains finite in both phases. Our results can be generally understood using 
bosonization and renormalization group results. The finite-size scaling of the charge Drude weight 
is well fitted by a polynomial function of the inverse system size in the metallic region. In the 
insulating region we find an exponential decay of the finite-size corrections with the system size 
and a universal relation between the charge gap Ac and the correlation length ^ which controls this 
exponential decay. 
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I. INTRODUCTION 

The transport properties of low-dimensional strongly 
correlated electron systems are currently a subject of 
great interest because of recent experimental observa- 
tions in quasi-onc-dimcnsional materials and because of 
their connection to the rapidly evolving field of nonequi- 
librium physics in strongly correlated quantum systemsii 
Most studies of transport properties in one-dimensional 
quantum many-body systems have been within the linear 
response theoryi^ One fundamental quantity is the Drude 
weight defined as the zero-frequency contribution to the 
real part of the conductivity. Thus a finite Drude weight 
implies ballistic transport, and it has been proposed as a 
criterion for distinguishing metallic and insulating phases 
in a Mott transition."^ 

Experiments on quasi-one-dimcnsional organic 
conductors^ show large deviations from the predictions 
of band theory. According to previous studies, the 
inter-site Coulomb repulsion plays a crucial role in these 
materialsi^ Therefore, the most simple effective model 
for their electronic properties is a one-dimensional ex- 
tended Hubbard model. In such onc-dimcnsional models 
of interacting fermions, the quasiparticlc concept breaks 
down, and the properties of the system do not resemble 
those of a Fermi liquid. Instead, low-energy excitations 
are made of independent elementary excitations for 
spins (spinons) and charge (holons)i^ii Moreover, the 
space- and time-dependent correlation functions display 
unusual power-law decays. Their exponents are not 
universal but depend on the strength of the interaction. 
Onc-dimensional metallic systems belong to the generic 



class of Tomonaga-Luttinger liquids (TLL)4i^ '^°i^^ 
Their characteristic quantities are the so-called TLL 
parameters Vp, v^, Kp, and Ka-. Here Vp and are the 
velocity of charge and spin excitations, respectively, and 
Kp and determine the algebraic decay of correlation 
functions. 

Recently, the discovery of the colossal magnetic 
heat transport in spin ladder materials such as 
(Sr,Ca,La)i4Cu2404i, where the magnetic contribution 
to the total thermal conductivity exceeds the phonon 
contribution substantially, has sparked interest in trans- 
port properties of quasi-lD spin modelsi ^^i^^'^"^ Under- 
standing the transport properties of theoretical models 
is of great importance for the interpretation of transport 
or NMR measurements but it is still an open problem 
for quantum systems involving many coupled degrees of 
freedomii^ Therefore, the development of methods for 
computing transport properties such as the Drude weight 
in strongly correlated systems and the investigation of 
theses properties are much anticipated. 

In this paper we study the charge and spin Drude 
weight of the one-dimensional extended Hubbard model 
at quarter filling using the density-matrix renormal- 
ization group (DMRG) method with periodic bound- 
ary conditions. This model is known to be "non- 
intcgrable" for general values of the parameters'^ and 
thus not amenable to an exact calculation of the Drude 
weight contrary to the Hubbard modelJ^ Investigations 
based on the g-ology^^ bosonization^ii^i2£ and the renor- 
malization group2ii2^ have provided analytic insight, 
particularly in the weak coupling regime. Both ex- 
act diagonalization22i2ii^i2Sii^ calculations and quantum 
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Monte Carlo simulations^^ have clarified a number of 
questions at intermediate and strong coupling. In the last 
decade the DMRG method has been successfully used to 
investigate many properties of one-dimensional strongly 
correlated lattice model a^^i'^° but a precise calculation of 
charge and spin Drude weights has not been reported yet. 
Therefore, we still lack a comprehensive picture of ballis- 
tic transport in the one-dimensional extended Hubbard 
model. 

Our paper is organized as follows. In Sec[TTl wc sum- 
marize some properties of the extended Hubbard model 
and its ground state phase diagram and introduce the 
charge and spin Drude weight. In Sec lIIH we explain the 
DMRG based method for calculating the Drude weight. 
Our results for the thermodynamic limit are presented in 
Sec HVH and the finit c-sizc scaling is discussed in detail 
in Sec lIVB] Finally, wc summarize our work in the last 
section. 



II. MODEL AND DRUDE WEIGHT 

Wc study the one-dimensional extended Hubbard 
model which is defined by the Hamiltonian 

H = Ht + Hu (1) 

/.cr 

Hu = U^ni^^ni^i + V^nini+i (3) 
I I 

where cj ^ (c/.<t) is the creation (annihilation) operator 
for an electron with spin a (=Ti i) site I = 1, . . . , L, 
'T-i.CT = c\ ^ci^a is the density operator, and ni = 7i;.|+?i/.j^. 
We use periodic boundary conditions throughout, t > 
is the nearest-neighbor hopping integral along the chain, 
[/ > is the onsite Coulomb interaction, and > is the 
nearest-neighbor Coulomb interaction. A quarter-filled 
band corresponds to a system with N = L/A electrons 
of each spin, a Fermi wavevector fcp = ^, and a Fermi 

velocity v-p = 2tsin(fcF) = ^f^t in the Fermi gas {U = 

V = 0). 

This model has been studied extensively by a variety 
of techniques. It is known to be "non-integrable" for 
general values of the parameters'^ on the basis of energy 
level statistics although exact results can be obtained in 
three limits (V = 0, t/ = +oo, V = +oo). For V = 0, 
the model becomes the regular Hubbard model. At 
quarter-filling, it is known to be mctallic'^^ with dominant 
2A:F-spin-density-wave (SDW) fiuctuations, and its low- 
energy excitations are of the TLL type. For U = -foo, 
the quarter-filled electron model is equivalent to a half- 
filled spinless fermion model which upon increasing V 
from zero undergoes a phase transition from a TLL 
phase to a 4A:F-charge-density-wave (CDW) insulator at 

V ~ 2t^ For V ~ +00, onsite electron pairs cannot 
move, while the unpaired electrons have the same ki- 
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FIG. 1: The phase diagram of the t-U-V model at quarter- 
filling determined from DMRG calculations in Ref. |30. In 
region (I), 1/3 < Kp < 1 and 2fcF-SDW correlations are dom- 
inant. Region (II) is characterized by dominant 4fcF-CDW 
correlation and 1/4 < Kp < 1/3. In region (III), triplet pair- 
ing correlations dominate because of Kp > 1. The region 
(I), (II), and (III) are metallic states while region (IV) is an 
insulating state with Akp charge ordering. 



netic energy as spinless fcrmions interacting with an in- 
finitely strong nearest-neighbor repulsion. This system 
has a Bethc ansatz solution^iiS^i^^iM. It is a 4A:f-CDW 
insulator for U > Uc = At whereas it is phase separated 
for U < Uc- In the weak-coupling limit {U,V ^ t), 
the model can be mapped onto a 5-ology model and in- 
vestigated using bosonization and renormalization group 
techniques^^iia^^^ 

The ground-state phase diagram of the quarter-filled 
extended Hubbard model for repulsive interactions was 
first determined using exact diagonalizationsi^ Recently, 
the precise ground-state phase diagram and the TLL 
exponent Kp have been obtained for a wide region of 
the C/y-parameter space using the DMRG method]^ 
These results are summarized in Fig. [TJ where four dif- 
ferent phases are represented: (I) A metallic phase with 
1/3 < Kp < 1 where the system has dominant 2A:f-SDW 
fluctuations, (II) a metallic phase with 1/4 < Kp < 1/3 
where the system has dominant 2fcF-CDW fluctuations, 
(III) a metallic phase {Kp > 1) where the system has 
dominant superconducting fluctuations, and (IV) an in- 
sulating phase (i.e. with a finite charge gap) where 
the system has a long-range ordered 4fcF-CDW. All four 
phases have gapless spin excitations j^^'^^i^^ Finite spin 
gaps have been reported in previous exact diagonaliza- 
tion studies for large V in phase (HI) but our DMRG 
calculations indicate that the spin gap vanishes in the 
thermodynamic limit at least for all V < lOt. On the 
metal-CDW transition line Kp = 1/4. It has been re- 
ported that this transition is of the Kosterlitz-Thouless 
type^^ and higher order scattering processes, including 
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the Sfcp-Umklapp scattering via the upper band around 
zbS/cp, play a crucial role in the weak-coupling theory of 
the phase diagram 7'^^i^°'^^ 



Let us consider external fields 
the kinetic energy 



and (pa which modify 



,i(0p+(T</)„)/L t 



h.c. 



(4) 



La 



(j)p is the magnetic flux threading the system and (f>a- is the 
magnetic flux given by a fictitious spin-dependent vector 
potentiali^ Time-dependent fields 0^ generate a charge 
{ji = p) or spin (/x = a) current (J^) in the system. 
Within the linear response theory the charge and spin 
conductivities have the forroi^i^i^i^i^ 



(5) 



where crj^'^s^^j) jg assumed to be regular at = 0. The 
coefficient Dp (Da-) of the 5- function S{uj) is called charge 
(spin) Drude weight for the ballistic transport, and is 
given by 



Sk + Sn 



1 

'2L 



is 



En — Eq 



(6) 



where tpn denotes an eigenstate of H with energy En 
and the ground state corresponds to n = 0. The first 
term (up to the prefactor —1/21/ ) is the total kinetic 
energy while the second term (up to a prefactor — 7r/2) is 
the total spectral weight of the incoherent part (tJ^°^([j) 
of the conductivity and describes the reduction of the 
Drude weight caused by incoherent scattering processes. 
The charge and spin current operators are Jp = j^j 
with local current operators defined from the continuity 
equation 



i [H,dpj] +jp,i+i -jp.i = 0, 



(7) 



where dp^i ~ ni and d^.i ~ ri;.! — n;,x. In our model, the 
precise form of j^^; is 

jp,i = ^Ui.aCi.a + li-c- (8) 

CT 

jaA = -itJ2'^ {4+l,a('l,a + ^■^■) ■ (9) 

(7 

The above spin current operator is different from the 
one used in spin modclsji^i^ Nevertheless, the spin 
Drude weights defined for an electron system or for a 
spin system have the same physical meaning because 
they characterize the response of the spin degrees of free- 
dom to the same external perturbation. The spin Drude 
weight Da is thus defined as a precise analog of the charge 
Drude weight Dp. Therefore, a value D^ > simply 
means that the system is an ideal spin conductor, so that 
the spin transport is not diffusivei^ 



Obviously, charge and spin Drude weight are equiva- 
lent in the non-interacting electron gas {U = V = 0) 
because the second term of © vanishes. Assuming that 
the low-energy excitation of our model can be expressed 
by the TLL theory, the charge and spin Drude weight 
can be represented byii 



Dp = UpKp/n 



(10) 



where Vp (va) is the renormalized charge (spin) velocity 
and Kp {K^) is the renormalized TLL exponent of charge 
(spin) mode. In our model, the renormalized value of 
is Ka = 1 because the system has a SU(2) spin-symmetry 
and there is no spin gap in the repulsive parameter region. 
Thus, irDa = Va- On the other hand, the behavior of Dp 
is more complicated because Vp and Kp depends on the 
interaction parameters. 



III. DMRG METHODS FOR DRUDE WEIGHTS 

A method for calculating the Drude weight ^ with 
DMRG was introduced several years ago^ but has been 
rarely used until now. In this section we briefly summa- 
rize our implementation of this numerical method and 
discuss some technical details. 

The first term of cq. (O can be easily calculated using 
the ground state DMRG method. The second-term can 
be calculated by targeting the correction vector \ip) which 
is solution of 



iH-Eo)\ij) = Jp\^Po). 



(11) 



The best implementation of this idea is a variational 
principle similar to the one used for the calculation of 
dynamical correlation fmictions4^i^ One considers the 
functional 



(12) 



If the ground state is not degenerate, this functional has 
a unique minimum for the quantum state which is the 
solution of (fTTI) . It is easy to show that the value of the 
minimum is related to the second term of eq. ^ : 



Wp{lprrdn) 



En — E{ 







(13) 



Our method consists in calculating the ground state and 
then minimizing this functional with DMRG. Note that 
this approach does not work for systems with a degener- 
ated ground state. Therefore, we always choose appro- 
priate system sizes L and numbers of electrons N to get 
a nondegenerate ground state. 

Another approach for obtaining the Drude weight with 
DMRG is to compute the dynamical current-current cor- 
relation function 



C,i,7, (w) = - (-^ol J,, 



1 



Eo-H 



— J^L iV'o) 

IT] 



(14) 
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using the dynamical DMRG (DDMRG) mcthod4^ii^ The 
imaginary-part of this quantity satisfies 



(15) 



which has been previously used to study the optical ab- 
sorption of various one-dimensional insulators including 
the extended Hubbard model at half filling4i The real 
part of the correlation function yields 



limReC^,,(0)=^ 



En — Eq 



(16) 



which can be expected from the Kramers-Kronig rela- 
tion with the f-sum rule of conductivity. Therefore, one 
can in principle calculate the second-term of ^ using 
DDMRG. However, a DMRG calculation using ^ and 
(fT3)) is faster and more accurate than a DDMRG calcu- 
lation of and (fT51) . In the first approach the error 
in the value of the minimum (V^min) is of the order of 

if we can calculate target states with an error of the 
order e < 1 within DMRG. With the DDMRG method 
the error in the real part of Cnniuj) is of the order of e 
(see the discussion in Refs. |42||43| ). 

As originally noted by Kohn;^ the Drude weight ^ 
can be calculated from the dependence of the ground 
state energy on the applied field using 



D„ 



L 



d^Eoi<f>^,) 



(17) 
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Therefore, one can calculate the Drude weight with 
DMRG using Eg. p?)) . This approach has been demon- 
strated on the spinless fermion model^^i^ but it requires 
treating complex Hamiltonians and performing a numer- 
ically delicate second derivative of the ground-state en- 
ergy with respect to 0^. Therefore, we have chosen the 
approach based on the variational principle p2|) and the 
equation p^ . 



FIG. 2: Contour map for (a) the charge Drude weight Dp and 
(b) the spin Drude weight Da in the i7y-parameter space of 
the extended Hubbard model at quarter filling. The bold 
line represents the boundary of the metal-insulator transition 



A. Thermodynamic limit 



IV. RESULTS 

We have carried out DMRG calculations for quarter- 
filled chains with periodic boundary conditions and 
lengths up to L = 60. The investigated system lengths 
are given by i = 8/ -I- 4 with integers Z > 1 so that 
the number of electrons of each spin {N ~ L/A) is odd. 
Thus the ground state has momentum P = and is not 
degenerate. We have kept up to m « 3200 density- 
matrix eigenstates in the DMRG procedure. The dis- 
carded weights are typically of the order 10~^ ~ 10^^ 
and the ground-state energy accuracy is ^ IQ^H. All 
energies and physical quantities are extrapolated to the 
limit m oo. 



We first discuss the DMRG results extrapolated to the 
thermodynamic limit {L — s- oo). The finite-size scaling 
is discussed in the next subsection. Contour maps of the 
charge and spin Drude weights are shown in Fig. [21 We 
can summarize our main results in five points, (i) Both 
Dp and Do- have their maximum 13^ = wf/tt = V2t/TT 
at the noninteracting point {U ^ V = 0) and decrease 
monotonically as a function of increasing U and V. (ii) 
Excepted for the noninteracting point, we observe a dif- 
ference between Dp and Do-- (This is due to the well- 
known spin-charge separation, which is a typical prop- 
erty in one-dimensional systems.) (iii) The charge Drude 
weight has no linear correction around the noninteract- 
ing point dDp/dU = dDp/dV = at U = V = 0. (iv) 
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FIG. 3: Charge Drude weight Dp and spin Drude weight 
as a function of U in the Hubbard model {V = 0). Lines show 
the exact results from the Bethe ansatz solution. 
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The charge Drude weight is larger than the spin Drude 
weight in the metallic phase, (v) Dp seems to be discon- 
tinuous at the metal-insulator transition. To understand 
these feature we will discuss the behavior of both Drude 
weights along the lines V = 0, U = 0, and U = IQt in 
more detail below. 

Our DMRG calculation results ioi V ~ can be 
compared with the exact Drude weights D™^'^' ob- 
tained from the Bethe ansatz solution of the Hub- 
bard modelii as shown in Fig. [31 Relative errors 

(•^cxact _ £)DMRG-) /£)oxact ^^^^^ f^^. ^^^^1 Sys- 

tem size L < 60 using up to 3000 density-matrix states. 
We note that the charge Drude weight is larger than the 
spin Drude weight for all [/ > in Fig. [3l The reduc- 
tion of both Drude weights for finite interactions can 
be understood qualitatively. The kinetic energy term 
Sk in eq.® is maximal for non-interacting electrons 
{U ~ V ^ 0) and decreases monotonically when U (or 
V) increases. The second term 5^ in eq.® equals at 
the noninteracting point and can only decrease to nega- 
tive values for ?7 > (or V > 0). (Note that S'p is not 
a monotonic function of U and V. It has a minimum at 
finite interactions as it converges to zero in the strong- 
coupling limit.) Therefore, the decrease of the Drude 
weight is due to both the suppression of the kinetic en- 
ergy and the appearance of scattering processes. 

In Fig. [3] we can see in both our DMRG results and 
the exact Bethe ansatz results that Dp has no linear cor- 
rection in U close to the noninteracting point {U = 0), 
whereas Da- seems to have a linear correction in U. This 
weak-coupling behavior has already been discussed li^i^ 
The renormalized value of is always K^- = 1 be- 
cause there is no spin gap for the parameters investigated 
here. Therefore, the lowest order correction to Do- yields 
Da = v-p/iT ~ U / {2tt^v-p), which is also consistent with 
the Bethe ansatz result 4^ The behaviour of Kp is more 
complicatedjii^ There is no first order correction from 



FIG. 4: (a) The charge Drude weight Dp and spin Drude 
weight Da as a function of V along the U = Q line, (b) The 
ratio Sp/S-K where and 5*^ are the first term and second 
term in eq.(|6]), respectively. Vertical lines show the boundary 
between the Luttinger liquid phases with dominant SDW (I) 
and pairing (III) fluctuations. Other lines are guides for the 
eyes. 



the interaction U and thus if we neglect the irrelevant in- 
commensurate 4fcF-Umklapp scattering^ and higher or- 
der corrections in Dp = VpKpj-n = v-p/ir. To un- 
derstand the decrease of Dp, the effect of the irrelevant 
4fcF-Umklapp scattering has to be taken into account. 
This correction is second or higher order in the interac- 
tion and explains the non-linear decrease of Dp. Note 
that, if we neglect the irrelevant 4fcF-Umklapp scattering 
but take higher order corrections in U into account (see 
Ref.12), Dp increases, which contradicts both our numer- 
ical results and the Bethe ansatz results. Thus, though 
those higher order corrections are important for quali- 
tatively understanding the phase diagram of this model, 
they are not sufficient for a quantitative analysis. 

In the strong coupling limit ?7 — > oo we expect the 
following behavior: Dp approaches the value t/yr because 
Vp 2tsin2fcF = 2t and Kp 1/2 while Da goes to 
because Va ~ 0{1/U). Therefore, along the ^ = 
line. Dp is larger than Da in both the U ^ and U —* 
oo limits and our DMRG results show that this relation 
holds also for all finite U > Q. 

Results for the U = line of the C/y-parameter space 
are shown in Fig. |D(a). Again we clearly see that both 
Drude weights decrease monotonically with increasing in- 
teraction and that the spin Drude weight is less than the 
charge Drude weight for all y > 0. In that case, however, 
it seems that both Drude weights have no linear correc- 
tion in V at the noninteracting point {U = V = 0). This 
is consistent with the weak-coupling theory which yields 
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FIG. 5: The charge Drude weight Dp and spin Drude weight 
Da as a function of V along the U = lOt line. The solid 
vertical line is the phase boundary between the metallic and 
insulating (IV) states determined by the calculation of Kp in 
Ref. [30I . The dashed vertical line marks the boundary be- 
tween the Luttinger liquid regions with dominant SDW (I) 
and CDW (II) fluctuations. 



no correction to VpKp or in first order in Vr^ The 
faster reduction of D^r with increasing V is due to the 
much stronger incoherent scattering processes for spin 
excitations than for charge excitations, which is demon- 
strated by the larger second term of eq. ([6]) as shown 
in Fig.lUb). 

We also find that S^/Sk increases in the Luttinger 
liquid phase (III) with dominant superconducting cor- 
relations, which occurs for V > 7.5t. This is an un- 
usual behavior for this "superconducting" phase because 
Sp/Sk = for a true superconductor 4^ In the V ^ 00 
limit, both Drude weights go to zero because the system 
becomes phase-separated at I/^ = 00 for U ~ Q. The 
lattice is decomposed in finite-size domains of singly oc- 
cupied or empty sites with an average electronic density 
< 1/2. These domains are separated by impenetrable 
immobile walls (the doubly occupied sites). There is a 
finite density of such pairs for small U because they re- 
duce the average density on the other sites below 1/2 
and thus allow them to gain kinetic energy. The charge 
Drude weight is zero despite the finite kinetic energy [i.e., 
S'k 7^ in ^] because charge motion is confined to finite 
domains by the infinite walls»^ Therefore, there are only 
incoherent contributions to the charge conductivity and 
Sp = —Sk- The spin Drude weight must also vanish for 
this reason and also because the infinite nearest-neighbor 
interaction V prohibits the formation of any nearest- 
neighbor electron pairs and thus the effective magnetic 
interaction is zeroj ^^'^^ 

Exact diagonalization studies^! have shown that the 
compressibility Kp increases with V in this "supercon- 
ducting" region and seems to diverge when approaching 
the phase separation regime for V 00. As within the 
TLL approach the compressibility can be characterized 




-log{Uf/t') 



FIG. 6: Spin Drude weight D^ vs - log(f/l/^/t^) in the strong 
coupling region. Solid lines are guides for the eyes. 



by Kp = 2Kp/TrVp, the divergence of Kp has been in- 
terpreted as a divergence of Kp. However, our DMRG 
calculations show that Dp = VpRp/ir decreases as V be- 
comes very large. Thus Vp goes to zero faster than Kp 
diverges toward 00 (if it diverges) and we conclude that 
the divergence of Kp is mostly due to the vanishing of 
the charge velocity for y ^ 00. Unfortunately, it is dif- 
ficult to compute Kp using the DMRG method close to 
the phase separation regime. Thus we cannot determine 
whether Kp diverges to 00 or converges to a finite value. 
Our results for the U = lOt line are shown in Fig.O As 

V increases, there is a phase transition from the metallic 
state to the 4A:f-CDW insulating state at Vc/t « 2.743 
In the metallic phase, both Drude weights decrease when 

V increases and we always find Dp > D„ as discussed 
previously. The spin Drude weight changes smoothly 
and remains finite for all values of the interaction pa- 
rameters that we have investigated. As seen in Fig. [5] 
the charge Drude weight approaches a finite value when 

V ^ Vc from the metallic side. As Dp = in the in- 
sulating phase for V > Vc, we conclude that the charge 
Drude weight jumps from a finite value to zero at the 
metal-insulator transition. Therefore, Dp = is a dis- 
continuous function of the interaction parameters. [Note 
that the finite-system Dp{L) calculated with DMRG are 
smooth functions of the parameters and from the finite- 
size scaling analysis alone one cannot determine whether 
Dp has a jump a,t V = Vc, see the next subsection.] The 
transition from the metallic to the insulating phase is 
believed to be of the Kosterlitz-Thouless type and to be 
caused by the S/sp-Umklapp scatteringj ^^'^° This scatter- 
ing is different from the usual 4A:F-Umklapp scattering 
which reduces the charge velocity Vp in renormalization 
group calculations and is always irrelevant . "^^'"^^ There- 
fore, both Vp and Kp are renormalized to a finite value 
when approaching the phase boundary from the metallic 
side in the renormalization group analysis. 
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Unfortunately, our simulations show that it is difficult 
to determined the phase boundary using the discontinu- 
ity in Dp. We believe that the direct calculation of Kp 
from correlation functions^ is a more efficient approach 
for determining phase boundaries within DMRG compu- 
tations because the extrapolation to infinite system size 
is less difficult than for the Drude weights. This will be 
discussed in the following subsection. 

In the strong-coupling limit of the insulating CDW 
phase (IV) electrons are localized and there is an ef- 
fective Heisenberg-like interaction between their spins. 
The effective exchange interaction J^ff between nearest- 
neighbor spins can be derived from perturbation theory 
and we obtain up to fourth order Jeft oc t'^/UV'^. In a 
Heisenberg model Du is proportional to the exchange in- 
teraction. Therefore we expect = Ct'^ /UV'^ in the 
strong-coupling limit of our electronic model, where C is 
an unknown constant. In Fig. [6] we show our numerical 
results for the strong-coupling limit on a double logarith- 
mic scale. One can clearly see a linear asymptotic behav- 
ior \og[D„/t) = log(C) - log(C/yVi3) for large U/t and 
V/t in agreement with the strong-coupling analysis. 



B. Finite-size scaling 

With DMRG we have been able to investigate much 
larger system sizes than in studies based on the ex- 
act diagonalization methodi^^i^ii^iS^iSi Finite-size effects 
are thus smaller than in these previous studies and in 
most cases we can perform a reliable finite-size scaling 
analysis of DMRG results for finite chains with periodic 
boundaries and extrapolated to the limit of infinitely long 
chains. Nevertheless, there are some difficulties close to 
phase boundaries as usual and it is always informative to 
examine finite-size corrections. 

Figure [7Ja) shows the spin Drude weight D„{L) as a 
function of inverse system length 1/ L for several param- 
eter sets. We note that Da{L) decreases monotonically 
and smoothly with increasing L and also that it converges 
to a finite value and seems to be a convex function of 
1/L for L ~* oo both in the metallic and the insulating 
phase. This finite-size scaling reflects the absence of gap 
in the spin sector and the complete separation of low- 
energy spin and charge excitations in the quarter-filled 
extended Hubbard model for U,V > 0. We can extrap- 
olate the Drude weight to the thermodynamic limit sys- 
tematically by performing a polynomial fit in 1/L. We 
find that the extrapolation is always very well behaved 
as seen from Fig. [Tj^a) . Especially, the relative errors are 
small, i.e. the difference between the extrapolated value 
Da and the value of Da{L) for the largest system size L 
computed with DMRG is small compared to D^. 

Figure. [7Kb) shows the charge Drude weight Dp{L) as 
a function of inverse system length 1/L for several pa- 
rameter sets in the metallic phase. As expected Dp{L) 
converges to a finite value for L ^ oo. As for the spin 
Drude weight we find that finite-size corrections are pos- 
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FIG. 7: (a) Spin Drude weight Do- as a function of inverse 
system length 1/L for several parameter sets in the metallic 
phase (filled symbols) and in the insulating phase (open sym- 
bols). Lines are polynomial fits, (b) Charge Drude weight 
Dp for several parameter sets in the metallic phase. Lines are 
polynomial fits. 



itive and depend smoothly on L and that Dp{L) seems 
to be a convex function of 1/L for L — > oo. Again such a 
finite-size scaling reflects the absence of gap in the charge 
sector in this part of the phase diagram. We can extrap- 
olate Dp{L) to the thermodynamic limit systematically 
using a polynomial flt in 1/L. This extrapolation is very 
accurate as demonstrated by the comparison with the ex- 
act results obtained from the Bethc Ansatz for the Hub- 
bard model (see Fig. 3). We note, however, that leading 
flnite-size corrections are linear in 1 / L in the U > V > 
region of the parameter space whereas they are of the 
order of 1/L^ in the Hubbard model.— In the V > U 
region we flnd that flnite-size corrections to Dp are much 
smaller than for U > V region and that the leading term 
seems to be again of the order L^^ as can be seen for 
U = 0,V = 5t in Fig. Ulh). 

In Fig. [H^a) we show the charge Drude weight Dp{L) 
in the insulating phase as a function of the system length 
L for several parameter sets. The linear behaviour ob- 
served in this semilog plot indicates that Dp{L) decreases 
exponentially with increasing L as in the Hubbard model 
at half flUingjii The correlation length ^ which char- 
acterizes this exponential decay, \og[Dp{L)/t] = — i/^, 
depends upon U and V. (We set the lattice constant 
a = 1.) Near the metal-insulator phase boundary one 
expects that should vanishes as the charge gap Ac. 
In particular, in the half-flUed Hubbard model it has been 



8 




6 

2 
0, 



1 ' 1 
-(b) 


1 


- • Ult=l 




. • Ult=% 




U/t=9 / 




" • U/t=\0 .^^^^^ 




.■•••"^1 1 1 


1 







0.1 



0.2 



0.3 



FIG. 8: (a) Semilog plot for the charge Drude weight Dp 
at U jt = 10 as a function of the system size L. Lines are 
exponential fits, (b) Charge gap Ac vs. inverse correlation 
length for \J jt = 10,9,8,7 and varying V. Solid lines 
are guides to the eyes. The dashed line is obtained from the 
Bethe ansatz solution^ in the limit I] jt = oo a&V approaches 
2t from above. 



shown that = Ac/2i = Ac/vf in the weak coupling 
limit U ^ t)^ We have determined the correlation length 
^ in our model from the slopes of log[Dp(L)/t] versus L. 
The charge gap is obtained as the difference of ground- 
state energies extrapolated to the thermodynamic limit 



Ac = lim Ac(L), 



(18) 
2£;l(0,0)]/2. 



Here El{N^ , N^) denotes the ground-state energy of a 
chain of length L with up-spin and down-spin 
electrons added or removed from the quarter-filled band, 
which can be easily computed using the DMRG method. 
Figure[5{b) shows Ac/t vs for several values of U and 
V. We see that varies linearly with Ac for small gaps 
while deviations are apparent when the gap is large. The 
product ^Ac tends to a universal value ^Ac w 13. St w 
9.76vp for Ac independently of U and V. This 
agrees with the relation between ^ and Ac derived from 
Bethe ansatz results^ for the effective spinless fcrmion 
model which corresponds to our model in the limit U/t = 
oo as F approaches 2t from above (see the dashed line in 

fig-iab)). 

It is also necessary to discuss the finite-size scaling of 
the charge Drude weight in the vicinity of the metal- 
insulator transition in more detail. In Fig. [5] we plot the 
charge Drude weight Dp{L) as a function of inverse sys- 
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FIG. 9: (a) The charge Drude weight Dp as a function of 
inverse system length 1/L at U/t = 10 for various V in the 
metallic phase (filled symbols, V <Vc ~ 2.74i) and in the in- 
sulating phase (open symbols, V > 14). Lines are polynomial 
fit. 



tem length 1/L for U = IQt and several values of V close 
to the critical value Vc ~ 2.74t determined in Ref . [sol. If 
we perform a polynomial fit in l/£, all results extrapolate 
to finite values for both phases while an exponential fit is 
meaningless. Thus from our numerical results alone we 
would conclude that Dp decreases smoothly as a function 
of V through the critical value Vc ~ 2.74t determined in 
Ref. [13 and vanishes only for V > V* , where V* is clearly 
larger than Vc- 

This failure of our finite-size analysis is easy to un- 
derstand. Obviously, if we assume a scaling behavior 
Dp{L) w ylexp(— L/^) in the critical region « Vc of 
the insulating phase and ^ diverges as V Vc, we need 
to treat exponentially large systems to observe the cor- 
rect finite-size scaling of Dp{L). Therefore, one has to 
use much larger sizes than the ones used here (up to 
L = 60) to determine the phase boundary from the ex- 
trapolated charge Drude weights. Unfortunately, it is 
very difficult to obtain accurate results for larger peri- 
odic systems using DMRG. We note that Dp{L) seems 
to be a convex function of 1 / L in the metallic phase but a 
concave function in the insulating phase. Although this 
finite-size behavior could be used as a criterion to dis- 
tinguish both phases in principle, it is not reliable gen- 
erally as it does not hold when V is kept fixed and U 
varies through the critical value Uc- Moreover, this is a 
transient behavior for intermediate values of 1 / L only as 
Dp{L) w j4cxp(— i/^) is also a convex function of 1/L 
for large enough i ^ ^ in the insulating phase. 

The Luttinger parameter Kp also extrapolates to a fi- 
nite value in the insulating phase if the available system 
sizes are too small;^ Nevertheless, the direct calculation 
of Kp from the finite-size scaling of correlation functions 
allows one to determine the metal-insulator phase bound- 
ary very accurately because this approach can be applied 
for open boundary conditions and thus with DMRG one 
can simulate systems which are one to two orders of mag- 
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nitude larger than with periodic boundary conditions. 

V. SUMMARY 

We have studied the transport properties of the t-U-V 
extended Hubbard model at quarter-filling by using the 
DMRG technique combined with a variational principle 
for the Drude weight. The contour map of the charge and 
spin Drude weights has been determined in the parame- 
ter space of the on-site Coulomb repulsion U and nearest- 
neighbor Coulomb repulsion V. We have found that (i) 
both Drude weights decrease monotonically with increas- 
ing Coulomb repulsion, (ii) the charge Drude weight is 
larger than the spin Drude weight in the metallic phase 
(Luttinger liquid), (iii) the charge Drude weight is dis- 
continuous across the Kosterlitz-Thouless transition from 
the metallic phase to the CDW insulating phase. 

We have also discussed the finite-size scaling and the 
extrapolation to the thermodynamic limit of our numer- 
ical data. In the insulating phase we find a universal re- 



lation between the charge gap and the correlation length 
which controls the exponential decay of finite-size cor- 
rections to the charge Drude weight. Unfortunately, we 
reach the conclusion that it is difficult to determine the 
phase boundary of a metal-insulator transition using the 
Drude weights calculated with DMRG because this ap- 
proach requires periodic boundary conditions which con- 
siderably reduce the performance of DMRG and thus the 
available system sizes for the finite-size scaling analysis. 
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